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ABSTRACT 


For  T  a  maximal  monotone  operator  on  a  Hilbert  space  H 
and  A  a  closed  subspace  of  H,  the  "partial  inverse"  TA  of 
T  with  respect  to  A  was  introduced.  TA  is  maximal  monotone. 

The  proximal  point  algorithm,  as  it  applies  to  TA,  was  shown  to 
result  in  a  simple  procedure,  the  "method  of  partial  inverses",  for 
solving  problems  in  which  the  object  is  to  find  x  e  A  and  y  e  AX 
such  that  y  e  T(x).  This  method  was  shown  to  specialize  to  give  new 
algorithms  for  solving  several  optimization  and  equilibrium  problems. 

An  algorithm  was  described  for  finding  a  feasible  point  for  a 
system  of  linear  inequalities.  If  the  solution  set  has  nonempty 
interior,  termination  was  shown  to  occur  after  a  finite  number  of 
iterations.  The  algorithm  is  a  projection  type  method,  .similar  to 
the  relaxation  methods  of  Agmon,  Motzkin,  and  Schoenberg.  It  differs 
from  previous  methods  in  that  it  solves  for  a  certain  "dual"  solution 
in  addition  to  a  primal  solution.  It  is  a  special  case  of  the 
method  of  partial  inverses. 

The  feasible  point  algorithm  was  shown,  for  inconsistent 
systems,  to  generate  a  sequence  converging  at  a  linear  rate  to  the  set 
of  least  square  solutions. 

A  primal-dual  decomposition  method  was  investigated  to  solve  the 
separable  convex  programming  problem.  Convergence  to  a  solution  and 
Lagrange  multiplier  vector  from  an  arbitrary  starting  point  was 
demonstrated.  Thi  method  was  shown  to  be  equivalent  to  the  method  of 
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partial  inverses.  In  the  nonseparable  case,  it  was  shown  to  specialize 
to  a  known  method,  the  proximal  method  of  multipliers.  Conditions  were 
provided  which  guarantee  linear  convergence  of  the  algorithm. 

For  families  of  nonlinear  programming  problems,  new  conditions  were 
established  which  guarantee  uniqueness  of  the  global  optimizer  to  be 
a  generic  property. 
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RESEARCH  SUMMARY 


Many  optimization  problems  are  equivalent  to  a  problem  of  finding 
a  zero  of  a  maximal  monotone  operator  T  on  a  Hilbert  space  H: 

(1)  to  find  x  e  H  such  that  0  e  T(x). 

The  known  "proximal  point  algorithm"  for  solving  (1)  takes  an 
arbitrary  starting  point  x^  e  H  and  a  sequence  (cn>  of  positive 
real  numbers  and  determines  a  sequence  (xn>  by  repeatedly  applying 
the  proximal  mapping: 

(2)  xn+1  -  (I*cnT)-l(xn). 

A  vast  array  of  problems  can  be  regarded  as  special  instances 
of  (1).  Examples  are  convex  or  linear  programming,  monotone 
complementarity  problems,  network  flow  problems,  variational 
inequalities,  and  systems  of  linear  equations  or  inequalities. 
Nevertheless,  the  known  uses  of  the  proximal  iteration  are  few.  By  far 
the  most  important  application  is  found  in  the  method  of  multipliers  of 
Hestenes  and  Powell.  This  algorithm  for  solving  convex  programming 
problems  was  known  for  some  time  before  Rockafellar  showed  it  to  be  a 
special  instance  of  the  proximal  point  algorithm.  Outside  of 
this,  there  are  few  serious  uses.  Although  it  is  theoretically 
possible  to  execute,  as  a  practical  matter  the  iteration  (2) 
can  be  performed  only  in  a  few  cases. 

There  are  many  problems  which  can  be  expressed  in  the  form  (1), 


are  not  easily  solved  by  the  proximal  point  algorithm,  but  have  the 
interesting  property  that  T  can  be  "decomposed"  into  simpler  monotone 
multifunctions  that  are  easy  to  proximate.  (To  "proximate"  means  to 
execute  the  rule  (2),  i.e.,  to  evaluate  the  proximal  mapping).  Our 
research  in  [1]  has  shown  that  new  algorithms  can  be  developed  which 
exploit  such  decomposition. 

The  principal  new  idea  that  enabled  us  to  extend  the 
powers  of  the  proximal  point  algorithm  and  derive  such 
decomposition  methods  is  our  notion  of  the  "partial  inverse" 
of  a  monotone  mapping.  If  A  is  a  closed  subspace  of  H 
and  B=AX ,  then  each  x  e  H  can  be  written  uniquely  as 
x  =  *A+xB  with  xA  e  A  and  x0  £  B.  If  T  is 
a  multifunction,  the  partial  inverse  of  T  with  respect  to  A 
is  the  multifunction  TA  :  H  *  H  defined  by  v  e  TA(u)  there 

exists  x  and  y  in  H  with  y  e  T(x),  u  =  *A+yB  and 
v  =  x„+y..  T.  is  (maximal)  monotone  iff  T  is  (maximal)  monotone. 

We  showed  that  TA  could  be  used  as  a  vehicle  for  introducing 
introducing  duality  into  and  solving  many  problems.  Often,  a  problem 
can  be  expressed  in  the  following  form  (for  suitable  choices  of  T 
and  A): 

(3)  to  find  x  e  A  such  that  there  exists  y  e  A 

with  0  e  TA< x+y ) , 

where  y  is  a  "dual  variable"  and  T  is  maximal  monotone.  If 
z  could  be  found  such  that  0  £  TA(z)  then  X=ZA  would  solve 


(3).  In  [1]  we  introduced  the  idea  of  solving  (3)  by  finding 
such  z  via  the  proximal  point  algorithm.  Theoretically,  this 
can  definitely  be  done;  the  fact  that  it  is  also  practical  was 
demonstrated  in  [1],  There  it  was  shown  that  this  leads  to  the 
following  iterative  procedure,  the  "method  of  partial  inverses"  to 
solve  (3): 


ALGORITHM  1  ("method  of  partial  inverses"  [1]). 
Initialize:  Choose  arbitrary  x^  e  A  and  y^  e  B. 
Iteration  k  (k=0,l,...): 

Proximal  step:  Find  x^,  y^  e  H  such  that 


xk+*k=xk 


Projection  step: 


and  1  £  T(<*kV  -  Uk)B)- 

ck  ck 

Let  =  (x^)A  and  yR+1  >  (y^B. 


In  [2],  we  applied  the  partial  inverse  approach  to  the  problem  of 
solving  a  system  of  linear  inequalities: 

(4)  to  find  x  e  R^  such  that  <x,u^>  <  b^,  i=l,...,n 

(where  0  /  u •  e  R^  and  b^  e  R).  Let  =  {x  :  <x,u^>  <  b-}  and 
C  =  Cx  n . . .  ncn.  It  is  a  straightforward  matter  to  apply  Algorithm  1 
to  solve  (4).  Our  principal  accomplishment  in  [2]  was  to  establish 
that  that  algorithm  actually  terminates  after  a  finite  number  of 


iterations,  provided  the  solution  set  has  nonempty  interior. 

The  algorithm  so  obtained  is  a  new  addition  to  the  family  of  known 


"projection"  methods.  These  solve  (4)  by  computing  a  sequence 
of  projections  onto  the  halfspaces  :=  {  x  :  <x,u^>  <  } . 

The  best  known  of  these  are  the  "relaxation"  methods  of  Agmon 
and  Motzkin  and  Schoenberg.  According  to  the  simplest  of  these,  a 
sequence  (x^)  is  generated  by  taking  x^^  to  be  the 
projection  of  x^  onto  the  furthest  halfspace  C^.  Another 
possibility  is  to  project  in  turn  in  some  fixed  cyclical  order 
onto  the  sets  C^.  In  another  variation  studied  by  Motzkin  and 
Schoenberg,  one  takes  x^+1  to  be  the  reflection  of  x^ 
across  the  furthest  hyperplane.  This  last  method  has  the  surprising 
property  that  a  solution  is  always  found  in  a  finite  number  of  steps 
if  the  solution  set  has  nonempty  interior.  In  a  more  recent  study  by 
Goffin,  classes  of  problems  (4)  were  identified  for  which  finite 
termination  in  the  relaxation  method  occurs  without  reflection. 
However,  the  reflection  method  is  the  only  one  of  these  known  methods 
for  which  finite  termination  always  occurs  under  the  sole  assumption 
that  the  feasible  set  has  nonempty  interior,  so  our  result  proving 
finite  convergence  of  the  algorithm  is  quite  surprising. 

The  algorithm  we  have  proposed  to  solve  (4)  is: 
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ALGORITHM  2. 

Start  :  Choose  arbitrary  Xq,  yg^ ,  ygn  e  R^ 

with  y01+. • -+y0n  =  0. 

Step  k  (k=0,l,...):  Compute 

xki  =  Pr°jc.(xk+yki)'  i=1'**"n 
yki  =  xk+yki‘xki  '  i=1 . n 


and  update 

x.  =  1  T?  .  x* 
k+1  —  ^ i=l  ki 


yk+l,i  yki  —  £j  =  l  ykj 
n  J  J 


i —1 f  •  •  • ,  n 


From  results  we  proved  in  [1]  it  follows  that  regardless  of  the  choice 
of  starting  values  xQ,  y^ , . . . , y0n, 


either  x^  -*>  x  and  yki 

-  yi 

with  x  e  C, 

yi 

normal  to  C^  at  x. 

and 

*1  *■■■+  m  °- 

or 

1  (xk+ykl' *  ‘  "xk+ykn)  1 

-+■  QO 

and  (4)  is  inconsistent 

Our  main 

result  from  [2]  regarding 

Algorithm  2  states  that 

termination  occurs  after  a  finite  number  of  iterations  if  the 
interior  of  C  is  nonempty.  More  precisely,  int(C)  /  0  implies 
for  some  k  that 


xk  =  xk+l  =  •  •  • 

0  =  ykl  =  yk+l,l  =  *  •  * 

0  =  i'kn  =  SVl.n  -  •  •  • 


In  [3],  we  investigated  the  behavior  of  Algorithm  2  in 
the  cases  where  int(C)  =  0  or  where  C  =  0.  If  int(C)  =  0  but 


C  t  0,  we  showed  that  the  sequence  (x^)  converges  to  a  solution  and 
that  the  distance  to  the  solution  set  approaches  zero  at  a  linear  rate. 
Even  in  the  case  where  C  =  0,  we  got  convergence  at  a  linear  rate  to 
the  set  of  least-square  solutions  of  (4),  and  (x^)  converges  to  one 
particular  least-square  solution.  More  precisely,  we  hav'  >roven  in  [3] 
the  following  (x  denotes  the  set  of  least  square  soluti'  ) 


THEOREM.  Let  sequences  (x^),  (y^),  (x£.),  ^yki^  -  produced 

by  the  Algorithm  2.  The  sequences  (x^)  and  (x£.)  (i  =  l,...,n) 
converge  to  limits 

1  •  xk  x» 

i 1  •  x^  -  x<;i,  i  =  l,...,n, 

where 

i.  x  =  1  Y-x'.  and  x'.  e  C- 

09  —  “1  ®  1  09 1  1 

n 

ii.  x;.  =  Pro jc^ ( xffl) 

x_  e  x  and  the  vector  (x' , -x  . . . . ,x'  -x  )  equals  the  vector  of 

09  /N  ©J,  09  ©n  09  * 

smallest  norm  in  the  set 

{ ( w, -  Yw; . w  -  Y  W ; )  •  W.£  C we  C  }  . 


This  vector  also  equals  the  element  v  of  smallest  norm  in  range(TA) 
For  the  sequences  (y^),  (y^),  one  ^as 

A 

iim  <*kl . ykn)/k  *  iim  <?kl . *kn}/k  “  -v- 

k-*-®  k>® 

In  [4],  we  applied  the  partial  inverse  approach  to  obtain  a 
new  decomposition  algorithm  for  the  solution  of  separable  convex 
programming  problems. 

The  convex  programming  problem  is 

(5)  to  minimize  fg(x)  subject  to  x  e  C  and 

f i (x)  <  0 ,  i  =  1, . . . ,m. 


We  assume  the  functions  f^  (i=0,...,m)  to  be  finite-valued 

convex:  Rd  -*•  R,  and  C  cz.  Rd  to  be  closed  convex.  In  the  separable 

case  of  ( 5) , 


(6) 


there  are  functions  f^  :  R  ^  ■*  R 

(d  =  d,+...+d  )  such  that  for  each  i, 
1  n 


i(x)  =  £"=1  f ij (xj ' 


(0<i<m,  l<j<n) 


di 

<x  “  <xi . V'  xj  e  R  > 

d  . 


and  there  are  closed  convex  sets  Cj  a  R  ^ 

such  that  C  =  C.X...XC  cr  Rd. 

l  n 


Our  method  is  closely  related  to  a  family  of  "dual"  methods  for 
the  solution  of  (5).  The  prototype  for  such  methods,  the  classical 
dual  approach,  involves  the  Laqranqian  function 


(xeC,  y>0 ) 


L(x,y)  =  fQ(x)  +  liml  y^f £<x) 

and  the  concave  dual  objective  function 

g(y)  *  inf  L(x,y)  (y>0). 

xeC 

The  dual  problem  to  (5)  is 

(7)  to  maximize  g(y)  subject  to  y>0. 

Under  mild  assumptions,  the  infimum  in  (5)  equals  the  supremum  in 
(7)  and  a  solution  y  >  0  to  (7)  exists.  Assuming  the  existence  of 
such  y,  it  is  known  that  x  solves  (5)  if,  and  only  if,  (x,y) 
is  a  saddle-point  for  L. 

In  a  typical  dual  approach,  one  minimizes  L(-  ,yk)  over  C  for 
a  sequence  of  values  yk  >  0,  obtaining  a  sequence  x^  which 
hopefully  converges  to  an  optimum  while  y^  converges  to  a  dual 
optimum.  Several  versions  of  this  strategy  have  been  suggested  in 
the  literature. 

One  very  valuable  characteristic  of  the  classical  approach  is  that 
it  leads  to  decomposition  algorithms.  For  the  separable  problem  (6), 

L(x,y)  =  fQ(x)  +  J.  ytf i(x) 


Thus  the  minimization  in  x  of  L(x,y)  over  C  decomposes  into  the 
separate  minimizations  of  Lj(xj,y)  =  ^0  j  ^ x j  ^  +  £i  ^i^ij^xj^  over  cj • 


u 


«q 
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This  "Lagrangian  decomposition"  has  been  exploited  by  numerous 


authors.  By  replacing  a  d-dimensional  constrained  problem  with  a 


sequence  of  less  constrained  ;  —oblems  of  lower  dimensions  d^,...,dn. 


it  offers  a  great  advantage.  Some  authors  have  investigated  ways  of 
applying  this  technique  in  situations  where  separability  is  absent, 
manufacturing  separability  by  the  device  of  replacing  functions 
with  their  linear  approximations  (linear  functions  are  always 


separable) . 


However,  the  classical  approach  has  several  potential  drawbacks. 


First,  for  some  y  >  0,  L(x,y)  may  fail  to  achieve  its  minimum  on  C. 
Worse,  g(y)  may  have  the  value  -<*>.  Such  values  of  y  must  be  ruled 
out,  so  the  dual  problem,  in  addition  to  the  nice  constraint  y  >  0,  has 
the  possibly  nasty  constraint  g(y)  >  -®.  It  is  possible  to  generate  a 
sequence  x^  failing  to  be  a  minimizing  sequence  for  (5)  even  with  yk 
being  a  maximizing  sequence  for  the  dual.  All  of  these  inconveniences 


can  be  ruled  out  by  imposing  appropriate  assumptions.  The  most 


serious  problem  is  that  one  is  severely  restricted  in  the  choice 


of  a  method  to  maximize  g.  Each  evaluation  of  the  function  g(y) 
requires  that  the  function  L(x,y)  be  minimized  in  x  over  C.  Thus 


any  method  requiring  many  evaluations  of  the  function  g  is 


impractical.  Fortunately,  the  minimization  in  x  of  L(*,y)  yields,  at 


no  extra  cost,  a  subgradient  for  g  at  y,  a  fact  which  motivates 


Uzawa's  method,  a  steepest  ascent  approach  to  the  maximization  of 


g.  The  Dantzig-Wolf e  algorithm  can  be  viewed  as  the  approach 


whereby  g  is  maximized  via  a  cutting  plane  method. 


One  attractive  strategy  for  avoiding  some  of  the  pitfalls 


%  A  A  L.V  A  M*  A  *"  /*  >•  A  .*  .  •  A  .  *  m’m  ■  A  .  '  A  A  »  *>  '»'  •  *,*  •  /  *  •  V  V 


“  •  *>  V  V  V  , 


SSI 


12 


of  the  classical  dual  approach  is  offered  by  Rockaf ellar ' s  proximal 
method  of  multipliers,  a  modification  of  the  multiplier  method  of 
Hestenes  and  Powell.  At  each  step,  one  has  x  and  y  at  hand 
and  minimizes  the  "augmented  Lagrangian" 

(8)  p(x)  =  f Q ( x )  +  1  | x-x | ^  +  E^max^fO,  f^(x)+y^} 

over  the  set  C  to  obtain  the  next  x.  The  strong  convexity  of  p(x) 
guarantees  the  existence  of  a  unique  minimum  and  the  next  multiplier 
y  is  chosen  according  to  a  simple  update  rule.  Global  convergence  is 
guaranteed  under  remarkably  weak  assumptions:  f^  finite  convex, 

C  closed  convex,  and  existence  of  a  solution-multiplier  pair. 

Unfortunatley ,  the  penalty  function  (8)  cannot  be  written  as 
a  sum  of  n  functions  p j  < x j ) ,  so  the  augmented  Lagrangian  approach 
does  not  directly  yield  a  decomposition  algorithm.  This  is  the 
principal  disadvantage  of  the  augmented  Lagrangian  approach  used  for 
decomposition  purposes.  Several  strategies  for  dealing  with  this  problem 
have  been  discussed  in  the  literature.  One  way  around  this  difficulty 
is  to  either  replace  p(x)  by  a  linear  approximation  or  rely  on  a 
method  that  minimizes  p(x)  that  uses  linear  approximations,  such  as 
Frank-Wolfe. 

In  the  nonseparable  case  (n=l),  our  method  specializes  to  the 
proximal  method  of  multipliers.  In  the  separable  case,  (8)  is 
replaced  by  a  function  of  the  form  £pj,  where  each  p^.  is  a  function 
only  of  Xj  and  looking  very  much  like  (8).  The  user  is  free 


to  minimize  the  strongly  convex  function  over  the  set  Cj  (this 

is  the  only  constraint  in  the  subproblem  minimization)  by  any  method 
desired.  The  advantages  of  the  proximal  method  of  multipliers  are 
retained.  One  still  has  existence  and  uniqueness  of  a  minimum  in 
each  subproblem  minimization,  and  global  convergence  to  a  solution 
and  multiplier  under  the  same  minimal  assumptions,  even  when  the 
minimization  is  performed  only  approximately  according  to  the 
stopping  criterion  we  provide.  The  update  rule  for  the  multipliers 
is  just  as  simple  as  in  the  proximal  multiplier  method. 

There  are  two  basic  approaches  to  the  solution  of  the  separable 
problem  in  a  hierarchical  or  multi-level  fashion.  These  are 
resource-directive  and  price-directive  methods.  In  the  resource 
directive  approach,  one  iteratively  determines  values  u^j  such 
that  the  solutions  of  the  problems 

(9)  min  fj(Xj)  subject  to  x^  e  and  f^Cx^)  <  u^ 

(j=l,...,n)  solve  (5).  In  the  price-directive  approach,  one 
iteratively  determines  "prices"  y^, . . . ,ym  such  that  the  optimal 
solutions  of  the  problems 

(10)  min  fgj(Xj)  +  £y^f^(Xj)  subject  to  xjeCj 


(j*l,...,n),  also  solves  (5). 

Our  method  differs  from  existing  methods  in  that  it  iteratively 
minimizes  a  function  depending  both  on  prices  yir...,ym  and 
resources  The  prices  converge  to  values  such  that  the 

solutions  to  (9)  are  solutions  to  (5)  and  the  allocations 
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converge  to  values  such  that  the  solutions  to  (10)  solve  (5).  The 
method  is  in  this  sense  both  price-  and  resource-directive. 

di 

Let  n  >  0  be  an  arbitrary  constant.  If  we  define  p  ^  :  R  J  -►  R  by 

f  1  ,  .1  , 

I  c  / _ t  \  .  I  _ »  I  ^  .  r  _ ^  f  a  r  /  —  i  \  _ _  'l  /  •  xr  ^ 


*  5*Iimax  (if  qj  £ 


pj(qj)“ 


(otherwise) 


the  decomposition  algorithm  we  have  introduced  is 


ALGORITHM  3. 

Initialize:  Start  with  arbitrary  x=(xlf...,x  )  e  Rd,  y  e  Rm, 
and  u  e  Rmn  such  that  T.u..  *  0,  i=l,...,m. 

*•]  i] 

Iteration  k  (k=0,l,...): 

Minimization  step:  For  j=l,...,n,  find  the  unique  qj  to 

minimize  Pj(qj)  (subject  to  the  implicit  constraint  ql  e  Cj  ) . 

Update:  Let  u^  =  max{u^  f  —  (qj ) }  and  then 

+  +  1  +  1 
xj  -  qj »  uij  -  u^j  —  uifc'  ^i  ~  ^i  +  nn  ^k  uik* 


Our  principal  result  regarding  convergence  of  this  algorithm  is: 


§  / 
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THEOREM. 


Let  Algorithm  3  be  applied  to  the  convex  programming 


problem  (5).  The  algorithm  is  equivalent  to  the  proximal  point 
algorithm  in  the  sense  that  (x+,u++y+)  =  ( I+(n3F)A)_1(x,u+y) ,  and 
P  is  a  certain  lower  semicontinuous  convex  function.  Suppose  the 
minimization  step  is  performed  accurately  enough  so  that 


in  step  k:  for  each  j,  dist (0 , 3p j ( q j ) )  < 


<Iek  <  co) 


holds.  If  the  generated  sequence  of  iterates  (x,u+y)  is 
unbounded,  there  exists  no  solution-multiplier  pair.  If  the  sequence 
is  unbounded  and  the  Slater  condition  is  satisfied,  then  the 
convex  programming  problem  (5)  has  no  solution.  If  the  sequence 
is  bounded,  then  x,  y,  and  u  converge,  respectively,  to  a 
solution,  Lagrange  multiplier,  and  optimal  allocation  for  the  convex 
programming  problem.  Convergence  occurs  at  a  linear  rate  provided 
the  strong  second-order  optimality  conditions  are  satisfied  for  the 
problem  (5). 

In  [5],  we  discussed  the  parameterized  family 

(Qp)  to  minimize  fg(x,p)  subject  to  x  e  C  and 

f^Xjp)  <  0  for  i  e  I  =  {l,...,s} 

=  0  for  i  £  J  =  {s+l,...,m}. 

By  a  "generic"  property  of  the  family  (Qp)  is  meant  a  property 
which  holds  for  all  problems  (Qp)  except  possibly  for  values  of  p 
in  a  subset  of  Lebesgue  measure  zero. 

Rockafellar  and  Fujiwara  have  (independently)  given  examples  of 


constrained  families  (Qp)  having  the  property  that  for  almost  all 
p,  (Qp)  either  has  a  unique  global  minimum,  or  no  global  minimum 
at  all.  Using  the  transversality  theorem,  we  showed  that  their 
results  can  be  extended  to  more  general  classes  of  parameter izat ions. 
Fuj iwara  considered  the  family 

(Q  )  to  minimize  f(x)-x-u  subject  to  g(x)=b+v 

U  f  V 

(f  :  Rn  -*  R,  g  :  Rn  -*•  Rm)  and  showed,  assuming  sufficient 
differentiability  of  the  functions  f  and  g,  that  for  all  v, 

(Q  )  has  at  most  one  global  solution  for  almost  all  u.  This 

U  r  V 

implies,  by  Fubini's  theorem,  that  for  almost  all  (u,v),  (Q  )  has 

U  r  V 

at  most  one  globally  optimal  solution.  He  proved  a  similar  result 
for  problems  where  the  equality  constraint  is  replaced  with  an 
inequality.  Fuj iwara  obtained  his  result  by  applying  a  theorem  of 
Araujo  and  Mas-Colell. 

Rockafellar  obtained  a  similar  result  for  the  more  general  family 

(Qw  u  v)  to  minimize  fg(v,x)+w»x  over  all  x  satisfying 

f.(v,x)+u^  <  0  for  i=l,...,s 

=  0  for  i=s+l , . . . ,m. 

Only  the  linear  perturbation  w  of  the  objective  function  has  a 
real  role  in  his  proof  of  global  uniqueness;  like  Fujiwara,  he  holds, 
the  other  parameters  fixed,  shows  the  result  to  hold  for  almost 
all  w  and  then  invokes  Fubini's  Theorem  to  show  that  the  result 
holds  for  almost  all  values  of  the  parameters  combined.  The 
principal  tool  he  uses  to  obtain  almost  sure  uniqueness  is  the 
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fact  that  a  convex  function  is  differentiable  except  on  possibly 


a  set  of  measure  zero  (although  convexity  of  the  functions  f^  is 
not  assumed).  The  Mas-Colell  theorem  employed  by  Fuj iwara  is  proved 
by  quite  direct  and  elementary  means.  In  a  more  recent  paper, 

Fuj iwara  has  established  a  generic  global  uniqueness  property  for 
constrained  problems  under  a  compactness  assumption  on  the  constraint 
set.  This  result,  however,  deals  with  "generic"  properties  which  hold 

for  all  functions  f^  in  an  open  dense  set  with  respect  to  the  strong 

2 

Whitney  C  -topology. 

Our  result  on  the  generic  global  uniqueness  property  can  be 
regarded  as  a  generalization  of  the  Mas-Colell  theorem  to 
constrained  problems.  Our  use  of  the  transversality  theorem 
simplifies  the  proof  considerably,  though  at  the  expense  of 
slightly  stronger  differentiability  assumptions.  We  considered 
families  of  the  form 

(Q  )  to  minimize  fn(x,p,q)  subject  to 

f i (x,q)  <  0,  i*l, . . . ,s, 

*  0 ,  i*s+l , . . . ,m, 

where  the  assumptions  on  f^  are  as  before  and  the  parameters  vary 
over  some  open  set  in  Euclidean  space.  The  criterion  we  used 
to  establish  generic  global  uniqueness  was: 

(11)  for  all  X£,  and  all  q,  the  function 

p  +  fg(x^,p,q)  -  fg(x2,p,q)  is  of  rank  one  at  all  p 
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